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We present a new formulation of the jamming phase diagram for a class of glass-forming fluids 
consisting of spheres interacting via finite-ranged repulsions at temperature T, packing fraction 
or pressure p, and applied shear stress S. We argue that the natural choice of axes for the phase 
diagram are the dimensionless quantities T/pa^, pet' jt^ and S/p, where T is the temperature, p is 
the pressure, H is the stress, a is the sphere diameter, e is the interaction energy scale, and m is the 
sphere mass. We demonstrate that the phase diagram is universal at low pa^ / e] at low pressure, 
observables such as the relaxation time are insensitive to details of the interaction potential and 
collapse onto the values for hard spheres, provided the observables are non-dimensionalized by the 
pressure. We determine the shape of the jamming surface in the jamming phase diagram, organize 
previous results in relation to the jamming phase diagram, and discuss the significance of various 
limits. 



I. INTRODUCTION 

Disordered solids made of molecules (glasses), mi- 
croscale particles, droplets, or gas bubbles in fluid (col- 
loidal glasses, emulsions, and foams), or solid macro- 
scopic particles in air (granular materials) can be flu- 
idized by raising temperature or kinetic energy, decreas- 
ing density or packing fraction, or applying a mechan- 
ical load such as a shear stress. These phenomena can 
be distilled into a "jamming phase diagram" as a func- 
tion of temperature, density, and shear stress [1 . Such 
a diagram describes whether the system is "jammed" 
{i. e. whether the relaxation time of the system 
exceeds some fixed, long time scale) or "unjammed" 
(z. e. whether the relaxation time is shorter than that 
time scale). While each system has a different jamming 
phase diagram, a jammed region generally exists at suf- 
ficiently low temperatures T, inverse packing fractions 
or shear stresses H. 

Recently, Xu et al. [2 found that the equilibrium dy- 
namical behavior of soft repulsive spheres is simplified 
considerably if one considers appropriate dimensionless 
quantities related to the relaxation time, temperature, 
pressure, and shear stress. In particular, it is useful to 
express the relaxation time of the liquid r, in terms of the 
time scale ^Jm/p(J^ where m is the particle mass, a is the 
particle diameter, and p is the pressure. This time scale 
characterizes the time for a sphere to move a distance 
equal to its diameter when accelerated by a typical com- 
pressive force of pcr^. Thus, it provides an estimate of the 
duration of a particle rearrangement. This choice is par- 
ticularly convenient because if one defines the dynamic 
glass transition by the criterion r^fpajm = constant, 
then in the low pressure limit, the glass transition is con- 
trolled solely by the ratio of temperature to pressure, 
T/pcr^, and is independent of the potential. Here, the 
low pressure limit corresponds to pa^ /e <C 1, where e 
characterizes the repulsive interaction energy scale. This 
limit also corresponds to the hard-sphere limit, where 



e ^ oo. Thus, all soft repulsive spheres behave as hard 
spheres when pcj^ /e 0. 

In this paper, we extend the analysis of Xu et al. [2] 
to describe the behavior under steady-state shear. The 
rheological behavior is characterized by the relation be- 
tween the shear stress S, and the strain rate 7. Following 
the approach of Xu et al., we consider the dimensionless 
variables S/p and j^/m/po■ and show that the steady- 
state shear rheology of soft repulsive spheres reduces to 
that of hard spheres in the limit where pcr^ /e 0. 

In addition, we suggest a new formulation of the jam- 
ming phase diagram in terms of the dimensionless tem- 
perature, shear stress, and pressure, T/pa^, S/p, and 
pcr^ /e. Here, the jamming surface is characterized by the 
criterion r^Jpa jm = constant instead of the standard 
criterion r = constant. One advantage of this formu- 
lation is that the diagram is universal for soft repulsive 
spheres in the limit of low pcj^ /e. In addition, the jam- 
ming transition of packings of frictionless spheres. Point 
J [3j, which controls many of the properties of pack- 
ings [4 , lies at the origin of the diagram. 

II. MODEL AND METHODS 

We consider a class of models of frictionless spheres 
with finite-range, repulsive interactions. In particular, 
we consider bidisperse spheres of mass m, half with diam- 
eter <j and half with diameter 1.4<j. The spheres interact 
through pairwise additive interaction potentials V{rij): 

V.ir,) = [^['-W) f-^'^<---, (1) 

[0 for Tij > CFij, 

where cr^j = ((Ji+<jj)/2 is the separation at contact. Note 
that Vb(r) is the hard-sphere potential. 

We study the steady-state shear rheology of three- 
dimensional systems by conducting non-equilibrium 
molecular dynamics simulations at fixed temperature T 
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and shear strain rate 7. Here, the Boltzmann constant 
kB IS set to unity. We define the temperature by the ve- 
locity fiuctuations relative to an imposed uniform shear 
gradient and impose a fixed shear strain rate using Lees- 
Edwards boundary conditions; the system is sheared in 
the x-direction with the shear gradient in the ^-direction. 
We use periodic boundary conditions in the z-direction. 
For hard spheres, we use an event-driven algorithm [5l|6] 
at fixed packing fraction cj) and periodically rescale ve- 
locities to keep the temperature within 1% of the desired 
value. For soft spheres, we use a conventional molecu- 
lar dynamics algorithm that numerically integrates clas- 
sical equations of motion. We employ Gaussian con- 
straints [IHS] to fix the instantaneous temperature T and 
pressure p. In all cases, we measure the average steady- 
state shear stress S and the relaxation time r defined 
by Ar,(r) = a/V3, where Ar,{t) = V((r,(t) - r,(0))2) 
is the root-mean-squared displacement in the vorticity 
direction. 



III. RHEOLOGY COLLAPSE IN THE 
LOW-PRESSURE LIMIT 

In this section we show that the rheology of repulsive 
soft spheres described by the potentials in Eq. [T] reduces 
in the low-pressure limit to the rheology of hard spheres. 
As discussed in the introduction, it is important to in- 
troduce the appropriate dimensionless quantities for de- 
scribing the shear stress S and the strain rate 7. The 
arguments of Xu et al. 0, suggest that it is most con- 
venient to use the pressure to obtain the dimensionless 
quantities S/p and j^/m/po■. For a given potential, di- 
mensional analysis tells us that we may write the depen- 
dence of the shear stress on the three control parameters, 
T, p, and 7, as a dimensionless function / of three inde- 
pendent dimensionless control parameters: 
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FIG. 1: (a) Shear stress, S, vs strain rate, 7, in standard 
molecular dynamics units for a = 2, two different values of 
T/pa^, and four different values of pa^/e. The filled sym- 
bols represent T/pa^ — 0.03 and the open symbols represent 
T/pa^ =0.1. The different shapes and colors represent dif- 
ferent pressures: pink diamonds are pcr^/e = 10~^"^; blue tri- 



angles are pa /e — 10 ; red squares are pa /e = 10 



and 



orange pentagons are pa /e 
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(b) Same data, made 



dimensionless by the pressure. For comparison, data for hard 
spheres are also shown in black at T/pa^ = 0.03 (black filled 
circles), and at T/pa^ = 0.1 (black open circles). 



Note that Eq. |2] isolates the interaction energy scale 
e in only one of the three control parameters, the di- 
mensionless pressure pcr^ /e. The combinations S/p and 
'^^m/pcF are familiar to the granular materials commu- 
nity [T0H12]. The dimensionless shear stress Tj/p is a 
macroscopic dynamic friction coefficient, while the di- 
mensionless strain rate '^^/ra/pa^ is typically called the 
inertial number and is understood physically as follows. 
As the system is sheared at constant pressure, it repeat- 
edly dilates and contracts. The dimensionless strain rate 
describes how fast the system is sheared relative to the 
time it takes for pressure to drive a dilated configuration 
into a close-packed configuration; it is the ratio of the 
contraction time to the shear time. 



In the limit pa^ /e ^ 0, we expect 




In Fig.[T] we demonstrate that the rheology for spheres 
with harmonic {a = 2) interactions collapses in the low- 
pressure limit onto the rheology for hard spheres, as 
expected. In Fig. [l] (a), we show the shear stress in 
standard simulation units, shear stress Sa^/e as a func- 
tion of strain rate ^^JmG^e. Data are presented for 
spheres with harmonic repulsions (a = 2 in Eq. fl]) at 
four different pressures, pa^ je = 10~^, 10~^-^, 10~^, and 
10~^"^, and two different temperatures, T/pa^ = 0.03 
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FIG. 2: Approach to the hard sphere Umit. (a) Dimension- 
less shear stress, S/p, vs pressure, p, at fixed dimension- 
less temperature, T/p — 0.1, and dimensionless strain rate, 
il\fv — O-l? for three different exponents, a. (b) Dimen- 
sionless relaxation time, rp^^^, vs pressure, p, for the same 
parameters, T/p = 0.1, / ^Jp = 0.1, and three different val- 
ues of a. In each plot, the horizontal line represents the value 
for hard spheres. 



and T/pcr^ = 0.1. In Fig. [T] (b), we show that for each 
value of T/pcT^ we can collapse the data from the four 
different pressures onto the hard sphere results by divid- 
ing the shear stress by the pressure and multiplying the 
strain rate by the time scale ^Jm|pG. All the data at 
the lower dimensionless temperature, T /pa^ = 0.03, col- 
lapse onto a hard sphere curve with an apparent dynamic 
yield stress at low strain rates, while all the data at the 
higher dimensionless temperature, T/pcj^ = 0.1, collapse 
onto a hard sphere curve with a linear viscous response, 
Ti/p (X "^^m/pcr^ at low strain rates. At higher strain 
rates, the system shear thins: T^/p grows more slowly 
than linearly with increasing "^^m/pcr. Note that for 
0.03 < T/pcj^ < 0.1 there is a continuum of curves, so 
the rheology changes continuously and is not described 
by two distinct branches, at least when plotted in this 
way. 

As shown in Fig. [2] (a), the rheology approaches a well- 
defined limit as pa^ /e for several different potentials. 
At fixed T/pa^ and '^^mjpG^ T^/p approaches a limit- 



ing value as pcr'^ /e 0. While dimensional analysis does 
not require that this limit should be the same for differ- 
ent potentials. Fig. [5] (a) shows that S/p approaches the 
hard-sphere value for several different exponents a (see 
Eq.[T]) as pa^ /e 0. Thus, the function F in Eq.js] char- 
acterizes the rheology of hard spheres. All potentials that 
vanish at a well-defined distance, such as the potentials 
of Eq. [T] reduce to the hard-sphere potential in the limit 
of zero overlap, so systems described by such potentials 
should exhibit hard-sphere behavior as pa Ve ^ 0. Eq.[3] 
therefore represents a universal limit of the rheology. 

Note that the data collapse of Fig.jljb) is not limited to 
the rheology but should apply to any observable quan- 
tity made dimensionless by the pressure. For instance, 
as we show in fig. [5] (b), the dimensionless relaxation 
time T^Jpa jm also approaches its hard- sphere value as 
pcr^/e 0. 



IV. NONLINEAR SHEAR RHEOLOGY OF HARD 
SPHERES NEAR THE GLASS TRANSITION 

The results for soft spheres at low pressures show that 
it is important to understand the rheology and relax- 
ation time of hard spheres in order to understand the 
corresponding behavior of a soft-sphere system. We have 
performed non-equilibrium molecular dynamics simula- 
tions for hard spheres under shear, with results that are 
summarized in Fig.[3j In these plots, connected lines rep- 
resent fixed values of T/pcr^. The results at fixed T/pa^ 
shown in Fig. [3] were obtained by interpolating from sim- 
ulations conducted at fixed packing fraction. I n Fig.[3| ^a), 
we plot the dimensionless shear viscosity rf^cr /pm^ de- 
fined as the ratio of the dimensionless shear stress S/p to 
the dimensionless strain rate '^^m/pcr^ as a function of 
E/p. In Fig. js^b), we plot the relaxation time r^jpajm 
vs S/p. The two plots are very similar: at high tem- 
peratures, both the viscosity and the relaxation time ap- 
proach limiting values at low shear stress, while for low 
temperatures, the viscosity and relaxation time increase 
without bound on the time scales accessible to our sim- 
ulations. 

Finally, Fig. [sj^c) shows the behavior of the packing 
fraction. For a given value of S/p, (\) decreases with in- 
creasing T/pcj^^ as expected; the kinetic energy pushes 
particles further apart at higher temperatures. At high 
temperatures, (j) is nearly independent of stress up to rel- 
atively high stresses, but at lower T /pa^ ^ (j) becomes more 
and more sensitive to S/p because the system dilates {(j) 
decreases) with increasing stress. 

We have cut off the data in Fig. [3] at a dimensionless 
shear stress of S/p ^ 0.5 since there is a dramatic change 
in behavior there: r^pcr/m increases while rf^fajprn^ 
Tj/pr^ and (j) decrease sharply. These changes are associ- 
ated with the onset of layering of the spheres in the plane 
perpendicular to the shear gradient, as indicated by long- 



4 




7t 




0.3 




TP 



0.1 







(b) 




\ . . . .zr::::^ 













10" 



10" 



10" 



FIG. 4: Ratio of shear viscosity to relaxation time for hard 
spheres. As in Fig. [3] each connected Une represents a differ- 
ent value of T/pa^ . (a) Ratio of shear viscosity to relaxation 
time in standard hard-sphere units, r/a^ /Tr vs S/p. (b) Ra- 
tio of dimensionless shear viscosity, r]y^a/pm, to dimension- 
less relaxation time, ry^pa/m. As in Fig. |3| the values of 



T/pa are: 0.02 (blue down triangles), 0.03 (red down trian- 
gles), 0.04 (black circles), 0.05 (blue circles), 0.06 (red circles), 
0.07 (black squares), 0.08 (blue squares), 0.09 (red squares), 
0.1 (black up triangles), 0.15 (blue up triangles), 0.2 (red up 
triangles), 0.25 (black pentagons), and 0.3 (blue pentagons). 



FIG. 3: Hard sphere results interpolated at prescribed val- 
ues of the dimensionless temperature, T/p. (a) Dimension- 
less relaxation time, ry^pa/m, vs dimensionless shear stress, 
S/p. (b) Dimensionless shear viscosity, riy^a/pm, vs S/p. (c) 
Packing fraction, 0, vs S/p. Each connected line is a different 
value of T/pa^: 0.02 (blue down triangles), 0.03 (red down 
triangles), 0.04 (black circles), 0.05 (blue circles), 0.06 (red 
circles), 0.07 (black squares), 0.08 (blue squares), 0.09 (red 
squares), 0.1 (black up triangles), 0.15 (blue up triangles), 
0.2 (red up triangles), 0.25 (black pentagons), and 0.3 (blue 
pentagons) . 

range order in the pair distribution function (not shown). 
This layering facilitates the shearing of layers relative to 
each other but impedes the mobility of spheres within a 
layer, causing the shear viscosity to decouple from the 
relaxation time. Such layering has been demonstrated 
to be an artifact of thermostats like ours that assume 
a linear shear profile [131 [M]. Thermostats that do not 
assume a linear profile yield similar results below the lay- 



ering transition but do not form layers at high stresses; 
instead, they exhibit shear thickening, where viscosity in- 
creases with shear stress. At the layering transition, our 
system shear thickens at fixed packing fraction but not 
at fixed T/p. Since the layering is an artifact of the ther- 
mostat, we focus on the isotropic phase below S/p = 0.5. 

The viscosity and relaxation time are often used in- 
terchangeably to locate dynamic glass transitions. This 
seems reasonable since Figs. [3]^a-b) show that their be- 
havior is qualitatively similar. In Fig. |4]^a) we show the 
ratio of viscosity to relaxation time in standard hard- 
sphere units, rja^/Tr. Although the ratio is nearly in- 
dependent of stress, the value of the ratio varies signifi- 
cantly with T/pa^. Fig.[4|b) shows the ratio of r]^/a/pm 
to r^/pcr/m^ v/^Pi pressure units. The result depends 
only weakly on T/pcj^ over the range studied. Thus, 
while viscosity can be used as a proxy for relaxation time 
for repulsive spheres at a given fixed pressure, care must 
be taken in comparing 77 to r if the pressure is varying. 
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The near-collapse shown in Fig. [4jb) implies that up 
to relatively high shear stresses, the hard sphere fluid 
behaves like a Maxwell fluid with a relaxation time that 
decreases with shear stress but a modulus rj/r that is 
independent of shear stress and proportional to the pres- 
sure. 



V. DIMENSIONLESS FORMULATION OF 
JAMMING PHASE DIAGRAM 

The simplicity of the behavior of repulsive spheres as a 
function of the dimensionless variables T/pa^, S/p, and 
pcr^/e suggests that it would be useful to recast the jam- 
ming phase diagram in terms of these variables instead 
of T, S, and packing fraction (j) [1 . Both sets of vari- 
ables span the same space of variables, but the choice 
{T/pcr^, T^/p^pa^ / e} has a distinct advantage. The plane 
at pcT^ /e = 0, spanned by T/pa^ and S/p, defines the 
hard-sphere limit, which is universal for finite-ranged re- 
pulsions. This choice of axes highlights the universality 
at pcr^/e = 0. Another advantage of this choice of axes 
is that the zero-temperature jamming transition. Point 
J [3], lies at the origin. 

It is instructive to first consider two-dimensional planes 
of the jamming phase diagram. Fig. [s] (a) shows con- 
tours of equal dimensionless relaxation time ry^pa/m^ 
separated by half decades, in the equilibrium plane at 
S/p = spanned by {T / pa^ ^ pa^ / for a soft sphere sys- 
tem with a potential described by Eq.[T]with a = 2. The 
specific shape of these contours depends on the potential. 
However, the downwards slope of the contours of equal 
T^/pa/m in Fig. [sja) is more generic; it refiects the fact 
that at fixed T/pcr^^ the dimensionless relaxation time 
decreases with increasing pa^ /e. As the potential softens 
or the pressure increases at fixed T/pcj^ ^ the amount of 
overlap increases. As a result, the soft spheres behave 
as hard spheres with a smaller diameter and relax more 
rapidly [15]. 

Notice in Fig. [sja) that the contours become more 
closely spaced as T/pa'^ decreases. This refiects the 
fact that spheres at fixed pcj^ /e are fragile glass- formers; 
that is, \og{T^p(j/m) increases faster than linearly with 
pcF^ /T. While the T/pa^-values of the contours decrease 
with increasing pcr^ /e^ the relative spacing between con- 
tours remains similar, indicating that the fragility, de- 
fined in terms of the dimensionless variables r^pa/m^ 
T/pcj^^ and pcj^ /e does not change significantly with in- 
creasing pa^ /e. If we instead define the fragility as the 
shape of the Xogir^pcr /m) vs 1/^ curve, we find that 
the system becomes less fragile-that is, \og(T ^ pa / m) 
vs 1/0 increases less steeply-as the softness, pcr^ /e^ in- 
creases, consistent with the interpretation of results for 
colloids of varied softness [16 . 

Fig. [sj^b) shows contours of equal r^fpajm^ separated 
by half decades, in the plane at pcr^ / e = spanned by 



{T/pcr^, S/p}. As we noted earlier, the limit pa^/e 
corresponds to the hard-sphere limit for any of the finite- 
ranged repulsions studied here. Therefore, the contours 
of equal r^/pa/m are universal in this plane. We con- 
structed Fig. [s] (b) from hard-sphere data. 

Fig. [sjc) shows the full three-dimensional jamming 
phase diagram in {T/pa^, T^/p^pcr^ /e} space for the same 
soft-sphere system with a = 2 in Eq. [l] We restrict 
ourselves to values of T/pa^ < 0.2, T^/p < 0.15, and 
pa^ /e < 0.2 to avoid the artifactual layering due to shear 
discussed in Sec. IV. Since we must interpolate to find the 
level sets of dimensionless relaxation time, it would be 
computationally expensive to construct entire surfaces. 
Instead, we draw contours where the surfaces intersect 
four planes: the equilibrium plane at T>/p = 0, the hard 
sphere plane at pcr^/e^ a plane at pa^/e = 0.1, and a 
plane at pa^ /e = 0.2. 

The diagram in Fig. [sjc) differs from the standard 
jamming phase diagram in two ways. First, the axes 
are different {{T/pa^^T^/p^pa^ /e} instead of {T, E,0}). 
Second, the surfaces represent contours of equal dimen- 
sionless relaxation time, T^Jpa~|m^ instead of contours 
of equal r or Tye/rna^ . Here we have shown con- 
tours for three values of r^pa/m^ separated by decades. 
The jamming surface for our system is given by the 
smallest surface closest to the origin and corresponds to 
T^pcF /m = 10~^ ; this is where we can no longer reach 
equilibrium on the time scale of our simulations. 

As expected, the dimensionless relaxation time in- 
creases monotonically as any combination of the three 
control variables {T/pa^ ^J^/p^pa^ /e} is reduced. The 
contours become more closely spaced as T/pa^^ S/p, or 
pa^/e decrease, signaling the rapid increase of relaxation 
time as the system approaches jamming. 

The jamming surface intersects the unstressed plane 
(S/p = 0) along a curve defined by {T/pcF^)g = 
fg{pa^/e)^ where {T/pcF^)g is the dimensionless dynamic 
glass transition temperature. This temperature depends 
on pressure, and decreases with increasing pa^ /e. How- 
ever, recent results indicate that the relaxation time of 
soft spheres can be mapped onto the relaxation time for 
hard spheres by approximating the soft sphere potential 
by a hard sphere potential with a smaller effective di- 
ameter [15 . In particular, any approach to the jammed 
surface, such as the black arrows of Fig. |5|d), can be 
mapped onto the approach for hard spheres, the purple 
arrow of Fig. [s^d) along pa^ je = S/p = 0. This im- 
plies that the physics of the hard-sphere glass transition 
governs the glass transition of all finite-range repulsive 
spheres. 

At nonzero shear stress, the jamming surface defines 
what could be considered either a stress-dependent dy- 
namic glass transition temperature or a temperature- 
dependent dynamic yield stress Y^y. For the hard-sphere 
system at pcr^ /e = 0, there is a nonzero dynamic yield 
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FIG. 5: Jamming phase diagram for a — 2. (a) Jamming phase diagram in the equiUbrium plane at S/p = spanned by 
\T I^G^ ^pa^ 1 1\. (b) Jamming phase diagram in the universal plane at pcr^/e = spanned by {T/pcr^, S/p}. Panel (b) was 
constructed using hard spheres, which we showed to be equivalent to soft spheres at pa^ /e — 0. In (a) and (b), we show 
contours of equal dimensionless relaxation time spaced by half decades: r^/pajm — 10°'^ (yellow up triangles), r^/pajm — 10 
(red diamonds), ry^pa /m = lO"*^"^ (pink pentagons), T\JpG jm — 10^ (blue squares), T\JpG jm — 10^"^ (orange down triangles), 
and Ty/pajm — 10^ (black circles). The contours are constructed by interpolation, (c) Full three-dimensional jamming phase 
diagram spanned by {T /pa^ ^T^/p^pa^ /e^ for a — 2. Lines represent logarithmically spaced contours of equal dimensionless 
relaxation time (top to bottom): r^/pajm — 10 (red), r^/pajm — 10^ (blue), and r^pa/m = 10^ (black). We show contour 
lines along four planes cut through the diagram: the equilibrium plane at S/p = 0, the hard-sphere plane at p = 0, a plane at 
p = 0.1, and a plane at p = 0.2. The surfaces are guides to the eye. (d) Schematic illustration of various paths in the jamming 
phase diagram. See text for details. 



stress for T/pa^ < (T/pcr^)^^, where (T/pcr^)^^ = 
fg{pa^ /e = 0) marks the dynamic glass transition for un- 
stressed hard spheres. At all pressures, the contours of 
constant relaxation time r^pcj jm must be quadratic in 
S/j9 at small S/p. This is expected from symmetry; the 
magnitude of the glass transition temperature should not 
depend on the sign of the shear stress. Thus, if the shape 
of the dynamic yield stress curve for hard spheres is given 
by E,/p = fHs{iT/pcj^)"S _ (T/pa^)), then f{x) ~ ^ 



for small x. A similar shape has been observed for con- 
tours of equal viscosity for a model metallic glass-forming 
liquid \T7\ . 

Now consider the low temperature limit, T/pa^ 0. 
In this limit, the shear stress S/p appears to satu- 
rate at a nonzero value set by the zero-temperature dy- 
namic yield stress S^o whose value depends on pressure: 
^yo/p = fo{po'^/^)^ where fo{x) > 0. In the hard- 
sphere limit, where pa^/e 0, the low-temperature 
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system has a nonzero dynamic yield stress (Fig. |5|d)): 
/o(pcr^/e) = fo{pcr^/e = 0) = constant. Thus, at low 
values of pa^/e for soft spheres, the athermal dynamic 
yield stress E^o must scale with pressure: H^o ~ p at 
T/pa^ = 0. This is consistent with granular experiments 
and simulations that find a macroscopic dynamic friction 
coefficient in the limit of low strain rate [TOHIIJ [181 Ej • 
We may use the jamming phase diagram to charac- 
terize other dimensionless observables besides ry^pa/m^ 
such as the packing fraction (/), as functions of the 
dimensionless control parameters {T/pa^ ,T^/p,pa^ /e}. 
For example, in the ffuid region each value of 
{T/pa^, T^/p^pa^ /e} corresponds to a certain value of the 
packing fraction; this is the equation of state. We find 
that ^ increases approximately linearly with 
agreement with the scaling for static sphere packings at 
T/pa^ = J^/p = The packing fraction decreases ap- 
proximately linearly with T/pa^, consistent with the idea 
that T/pcr^ controls the amount of free volume. It also 
decreases nearly linearly with S/p, but with a smaller 
coefficient. 

In the jammed region below the jamming surface, 
many studies have focused on the behavior at T/p = 
Ti/p = as packing fraction is reduced toward Point J, 
the point at T = E = and = <pc ^ 0-64 where most 
static sphere packings lose mechanical stability [4 . This 
path corresponds to decreasing pa^/e toward the origin 
along T/pcr^ = E/p = 0, as shown by the red arrow in 
Fig. [5|d). Brito and Wyart j2QH22] have argued that a 
hard-sphere glass {i.e. a system at pa^/e = below the 
gray surface in Fig. jsjd)) can effectively be mapped onto 
a jammed, athermal soft-sphere one at time scales short 
compared to the time between rearrangements due to ag- 
ing. The effective coordination number and vibrational 
properties of the system along the path described by the 
green arrow at pa^/e = S/p = can be mapped onto 
those along the red arrow at T/pa^ = S/p = [2QH22] . 
It is likely that the short-time-scale properties of repul- 
sive spheres along any path that approaches the origin of 
the jamming phase diagram, such as the paths described 
by the white arrows in Fig. jsjd), are controlled by the 
same physics. 

However, in the jammed region the value of the pack- 
ing fraction at a given state point {T/pcr^ ^T^/p^pcr^ /e} 
is not given by an equation of state but becomes history 
dependent. For example, if T/pa^ is decreased slowly at 
T^/p = pa^ /e = 0, then the resulting value of (j) may be 
higher than if T/pa^ were decreased quickly to the same 
final value. The equation of state becomes multi-valued; 
at least one extra parameter, for instance the packing 
fraction, must be specified in order to determine the state 
of the system. Within the mean-field theory of random 
first-order models, glass states are uniquely determined 
by the jamming state diagram obtained by adding one 
more parameter to the jamming phase diagram |23j, but 
in general the state of the glass may depend in more de- 



tail on the history. 

Note that because the jammed state is history- 
dependent, the value of (j)c that corresponds to the jam- 
ming transition at the origin of the jamming phase 
diagram, {T/pa^ = O^^/p = 0,pcr^/e = 0}, is 
not unique. However, the approach to the origin in 
{T/pa^,T>/p,pa^ /e} space is well defined. This is an- 
other advantage of this formulation of the jamming phase 
diagram. 

Finally, we note that there is one region of parame- 
ter space that cannot be explored in the reformulated 
version. All systems at strictly zero temperature, zero 
applied stress, and packing fractions below the jam- 
ming transition, (j) < (j)c^ have zero pressure, so these 
states cannot be uniquely determined by the values of 
{T /pa'^ .,T^/p.,pcF^ / e}. However, as long as the tempera- 
ture or stress are not strictly zero, the pressure becomes 
nonzero, and the system can be represented on the refor- 
mulated diagram. 

VI. SUMMARY 

We have demonstrated that the jamming transition 
that occurs as a function of some combination of tem- 
perature, pressure or packing fraction, and applied me- 
chanical load is conveniently described in terms of the 
dimensionless quantities T/pa^, pcr^/e, and S/p. Such 
a formulation defines the three-dimensional state space 
spanned by T, p, and H as a product of the hard-sphere 
plane at pcj^ /e = and the unstressed plane at H/p = 0. 
One advantage of this formulation is that all repulsive 
spheres act as hard spheres near pcr^/e = 0, so the jam- 
ming surfaces for different repulsive spheres collapse on 
that plane. A second advantage is that any repulsive 
sphere system will undergo a dynamic jamming transi- 
tion as any combination of {T/pa^ ^T^/p^pa^ /e} is de- 
creased before reaching Point J at the origin. 

While the jamming surface is universal at low pcr^/e, 
it depends on the shape of the interaction potential at 
higher pa^/e. At zero shear stress, it has recently been 
shown that the relaxational dynamics of repulsive spheres 
can be mapped onto the dynamics of hard spheres using 
only structural information [TF, even for high values of 
pcr^/e. These results suggest that, at least for repulsive 
spheres, the entire {T / pa^ , pa^ / e} plane at S/p = can 
be mapped onto the hard-sphere path denoted by the 
purple arrow in Fig.jsj Previous work suggests that a rea- 
sonably well-defined effective temperature controls the 
behavior of soft spheres under shear stress [24H28] . One 
open question is whether the relaxation time of a system 
at nonzero S/p, corresponding to some value of Teff/pa^, 
can be mapped onto the equilibrium hard-sphere behav- 
ior at the corresponding value of T/pa^. 

Finally, we note that the results presented here apply 
to repulsive spheres only. In particular, the result that 
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the jamming phase diagram is universal at low pcj^ /e does 
not hold for systems with attractions. In order to de- 
scribe jamming of molecular liquids, we must also under- 
stand the effect of attractive interactions on the jamming 
phase diagram. 
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